TECHNICAL  REPORT 


THEORETICAL  PREDICTION  OF  THE  FLOW 
IN  THE  WAKE  OF  A  HELICOPTER  ROTOR 

PART  1  ■  DEVELOPMENT  OF  THEORY  AND 
RESULTS  OF  COMPUTATIONS 


By:  Peter  C 


- 


rimi 


CAL  No.  BB-1994-S-1 


Prepared  for: 

U.S.  Army 

Ballistic  Research  Laboratories 
Aberdeen  Proving  Ground,  Maryland  21005 


Final  Report  -  Part  1 

Contract  No.  DA30-069*AMC-645(R) 

September  1965  _ 

/ 


*  t 
t  , 


c/,Gc  0’7  S 


i  %  ':W' 


\\ 

jiitt'TGBy 

b 


CORNELL  AERONAUTICAL  LABORATORY,  INC. 


DISCLAIMER  NOTICE 


THIS  DOCUMENT  IS  BEST  QUALITY 
PRACTICABLE.  THE  COPY  FURNISHED 
TO  DTIC  CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 

I 


CORNELL  AERONAUTICAL  LABORATORY,  INC. 
BUFFALO,  NEW  YORK  14221 


THEORETICAL  PREDICTION  OF  THE  FLOW 
IN  THE  WAKE  OF  A  HELICOPTER  ROTOR 

PART  I 

DEVELOPMENT  OF  THEORY 
AND  RESULTS  OF  COMPUTATIONS 


By: 

PETER  CRIMI 


CAL  NO.  BB- I994-S- 1 

FINAL  REPORT  -  PART  i 
CONTRACT  NO.  DA30-069- AMC-645  (R) 

SEPTEMBER  1965 


Prepared  for: 


U.S.  ARMY 

BALLISTIC  RESEARCH  LABORATORIES 
ABERDEEN  PROVING  GROUND,  MARYLAND  21005 


SUMMARY 


A  study  was  performed  with  the  objective  of  predicting  the  time-varying  flow 
in  the  vicinity  of  a  helicopter  in  hovering  or  forward  flight.  An  analytical 
model  was  formulated  which  represents  the  rotor,  its  wake  and  the  fuselage 
of  the  aircraft.  The  model  for  the  wake  has  been  simplified  in  a  manner 
consistent  with  flow  visualization  studies  by  representing  the  wake  as  tip 
vortices  emanating  from  each  blade.  The  wake  vortices  are  assumed  to  be 
free  to  convect  under  the  influence  of  the  blades,  the  fuselage  and  the  wake 
itself.  The  fuselage  has  been  replaced  by  a  surface  distribution  of  sources, 
the  strength  per  unit  area  being  determined  by  assuming  that  the  fuselage  is 
subjected  to  a  uniform,  steady  free  stream. 

A  digital  computer  program  v/as  written  which  implements  the  model  developed. 
Flow  calculations  were  performed  which  agree  well  with  measurements  of 
the  time  average  of  the  flow  in  the  wake  of  a  two-bladed  rotor.  Extensive 
calculations  were  then  carried  out  which  define  the  flow  in  the  vicinity  of  a 
UH-1B  helicopter  for  a  range  of  flight  conditions  from  hover  to  high-speed 
cruise. 

Upon  analysis  of  the  results,  it  was  concluded  that  the  model  developed 
provides  a  valid  representation  of  the  unsteady  flow  in  the  vicinity  of  a 
helicopter.  It  was  pointed  out  that  the  formulation  is  somewhat  restrictive, 
in  that  the  flow  in  the  rotor  plane  and  immediately  adjacent  to  the  fuselage 
cannot  be  accurately  predicted.  However,  with  suitable  refinements 
incorporated,  these  regions  could  be  investigated  as  well. 
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FOREWORD 


The  work  reported  herein,  performed  between  September  1964  and  September 
1965,  was  accomplished  by  the  Cornell  Aeronautical  Laboratory,  Inc.  (CAL), 
Buffalo,  New  York  for  the  Director  of  Ballistic  Research  Laboratories,  (BRL) 
Aberdeen  Proving  Ground,  Maryland.  The  research  effort  was  performed 
under  Contract  DA  30-069-AMC-645(R)  and  was  monitored  for  BRL  by 
Mr.  Thomas  Coyle  as  Technical  Supervisor.  Dr.  Peter  Crimi  of  CAL 
conducted  the  study  and  received  assistance  from  Mr.  Alexander  Sowydra 
during  the  development  of  the  mathematical  model  and  Mr.  Harvey  Selib  for 
the  digital  computer  programming. 

This  document  is  Part  1  of  the  final  report  under  the  contract.  It  describes 
the  development  of  the  theory,  discusses  the  results  of  the  computation,  and 
provides  a  comprehensive  discussion  of  the  work  performed  under  the  contract. 
Part  2  of  the  final  report  describes  the  formulation  and  application  of  the  rotor 
wake-flow  computer  program  and  is  of  use  primarily  to  those  who  plan  to  use 
the  digital  computing  program. 

CAL  Report  Numbers  have  been  assigned  as  follows: 

BB-1994-S-1,  THEORETICAL  PREDICTION  OF  THE  FLOW  IN 
THE  WAKE  OF  A  HELICOPTER  ROTOR,  Part  1  -  Development 
of  Theory  and  Results  of  Computations 

BB-1994-S-2,  THEORETICAL  PREDICTION  OF  THE  FLOW  IN 
THE  WAKE  OF  A  HELICOPTER  ROTOR,  Part  2  -  Formulation 
and  Application  of  the  Rotor  Wake  Flow  Computer  Programs 
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SYMBOLS 


Oj  vortex  core  radius 

number  of  rotor  blades 

/*•  designation  for  iih  wake  reference  point 

R  rotor  radius 

r  position  vector  of  a  point  in  the  flow 

{V*  »  components  of  fluid  velocity  in  the  directions  of  x.  y  and 

respectively,  nondimensionalized  by  All/? 

(VK,  time  averages  of  V^,  Vy  and  Vj  respectively 

wake  and  blade  induced  contributions  to  and  V j.  respectively 

free-stream  velocity 

v'  induced  downwash 

v*  time  average  of  \r 

v'o  momentum  value  for  the  time  average  of  v 

W  aircraft  weight 

(x>>y*})  rectilinear  coordinates  with  origin  in  the  tip-path  plane, 
nondimensionalized  by  R 

&r  inclination  of  the  free  stream  to  the  tip -path  plane 
r  circulation 

A  blade  loading  parameter;  A  =  4 W/(ttzNb /OllzR¥) 

ju,  advance  ratio;  ju,  -  Vf  J [CiR) 

p  air  density 

cr  source  strength  per  unit  area 

f  rotor  azimuth  angle 

jfl  rotor  angular  speed 
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GLOSSARY  OF  TERMS 


Circulation: 

Convection: 

Roll-up  of  the  Wake; 

Self-induced  Velocity: 

Vortex: 


Vorticity : 


The  line  integral  about  a  closed  curve  of  the 
component  of  fluid  velocity  tangent  to  the  curve. 

The  process  of  transport  of  fluid  particles. 

The  concentration  of  wake  vorticity  into 
isolated  vortices. 

The  contribution  to  the  velocity  of  convection  of 
a  vortex  from  the  vortex  itself. 

A  concentration  of  rotational  fluid  (i.  e. ,  having 
nonzero  vorticity),  resulting  in  a  swirling  flow; 
bound  -  constrained  to  move  in  a  prescribed 
manner;  free  -  allowed  to  convect;  root  -  vortex 
emanating  from  the  blade  root;  sheet  -  a 
distribution  of  concentrated  vortices  in  the  form 
of  a  surface;  tip  -  vortex  emanating  from  the  blade 
tip. 

The  curl  of  the  fluid  velocity;  a  measure  of 
rotationality. 


Wake  of  a  Blade: 


The  vortex  elements  created  by  the  passage  of 
the  blade  through  the  fluid. 


1.  INTRODUCTION 


A  method  for  accurately  predicting  the  time -varying  flow  in  the  vicinity  of  a 
translating  or  hovering  helicopter  would  be  of  considerable  value  as  both  an 
engineering  tool  and  an  aid  to  research  and  development.  Possible  areas  of 
application  are  the  determination  of  harmonic  airloads  on  the  blades,  fuselage 
or  auxiliary  lifting  surfaces,  investigation  of  rotor  transient  response  and 
flutter,  vibration  and  fatigue  analyses,  and  study  of  downwash  impingement 
and  debris  entrainment  in  ground  effect.  The  method  of  flow  prediction 
which  is  reported  here  is  presently  limited  to  application  in  only  certain 
of  these  areas,  but  it  appears  that,  with  further  development,  the  scope  of  its 
applicability  might  be  considerably  broadened. 

Direct  analysis  of  the  flow  in  the  wake  of  a  rotor  for  arbitrary  flight  conditions 
is  a  formidable  task,  particularly  if  fuselage  effects  are  important  in  the 
region  of  interest.  The  dependent  variables  are  nonlinear  functions  of  three 
space  variables,  time  and  each  other,  with  no  means  available  for  linearization. 
Some  results  have  been  obtained  previously  by  retaining  the  complex  structure 
of  the  flow  but  specifying  the  wake  geometry.  For  example,  good  agreement 
with  experimental  measurements  of  blade  loading  were  obtained  by  calculating 
the  flow  at  the  rotor  plane  in  this  way  (Reference  1).  The  time  average  of  the 
induced  downwash  in  certain  limited  regions  was  fairly  well  approximated 
using  a  similar  assumption  as  to  wake  geometry  (Reference  2), 

The  objective  of  the  study  reported  here  is  to  predict  the  time -varying  flow 
at  an  arbitrary  point  in  the  wake  of  a  helicopter  rotor.  If  the  wake  geometry 
were  assumed  known,  unacceptably  large  errors  would  result.  Instead,  the 
wake  structure  has  been  simplified  in  a  manner  consistent  with  the  findings 
of  flow  visualization  studies  (see  Chapter  2).  The  wake  can  then  be  allowed 
to  convect  under  the  influence  of  the  blades,  the  fuselage  and  self-induced 
effects  by  constructing  a  numf  rical  analogue  with  the  aid  of  a  digital 
computer. 


In  developing  the  theory,  consideration  has  been  limited  to  aircraft  having  a 
single  rotor.  No  restriction  has  been  placed  on  the  number  of  blades  which 
the  rotor  may  have,  however.  The  aircraft  has  been  assumed  to  be  in  steady 
level  flight  or  hover,  out  of  ground  effect,  and  the  forward  speed,  rotor 
rotational  speed  and  loading  have  been  assumed  to  be  within  the  limits  of 
presently  operational  single-rotor  helicopters.  Also,  in  representing  the 
fuselage,  the  flow  has  been  assumed  to  be  attached  and  noncirculatory  and 
the  fuselage  geometry  has  been  somewhat  idealized  by  omitting  appendages 
such  as  landing  skids. 

The  results  of  the  study  are  reported  in  two  parts.  The  development  of  the 
theory,  the  basic  formulation,  the  major  assumptions  for  numerical  analysis 
and  the  results  of  computations  are  presented  here.  The  second  part,  contained 
in  a  separate  document,  concerns  the  detailed  information  related  to  the  coding  ' 
and  operation  of  the  digital  computer  programs  which  implement  the  theory. 
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2. 


DISCUSSION  OF  THE  PHYSICAL  FLOW 


It  is  desired  to  analytically  define  the  flow  in  the  vicinity  of  a  helicopter  in 
translational  or  hovering  flight  out  of  ground  effect.  In  order  that  a  rational 
framework  for  the  analytical  model  may  be  constructed  the  essential  features 
of  the  physical  flow  should  first  be  outlined. 

Generally,  there  are  three  primary  contributions  to  the  flow  at  a  given  point 
relative  to  the  aircraft.  Specifically,  the  rotor  blades,  their  wake  and  the 
fuselage  all  affect  the  fluid  velocity.  These  three  effects  are,  of  course, 
interrelated  in  an  extremely  complicated  and  nonlinear  manner,  but  for  the 
purposes  of  this  discussion  they  can  be  considered  separately. 

The  rotor  blades  may  be  regarded  as  wmgs  of  very  high  aspect  ratio  in  a 
harmonically  varying  stream.  While  the  wake  and  fuselage  do  alter  the 
loading  on  the  blades,  the  total  lift  at  a  given  azimuth  position  of  the  rotor 
must  still  be  approximately  equal  to  the  weight  of  the  aircraft.  This  implies 
that  the  circulation  about  a  given  blade  must  vary  approximately  sinusoidally 
with  azimuth.  The  spanwise  variation  of  circulation  about  a  blade  is  also  a 
function  of  wake  and  fuselage  effects,  and  depends  on  blade  twist  and  planform 
as  well.  Results  of  theoretical  computations  of  harmonic  airloads  (Reference 
1),  which  agree  well  with  experimentally  obtained  data,  indicate,  though,  that 
the  spanwise  variation  of  circulation  generally  does  not  differ  a  great  deal 
from  an  elliptic  distribution. 

The  wake  of  each  blade  is  generated  as  a  thin  sheet  of  vortical  fluid.  The 
vorticity  in  1  lis  sheet  has  components  in  the  spanwise  direction  (shed  vortices) 
and  in  the  chordwise  direction  (trailing  vortices).  It  has  been  observed  from 
flow  visualization  studies  that  this  sheet  rolls  up  within  a  few  chord  lengths  to 
form  a  single  concentrated  vortex  trailing  from  near  the  tip  of  the  blade  (see, 
for  example,  Reference  3).  While,  at  least  in  theory,  a  vortex  should  also 
trail  from  the  vicinity  of  the  root,  it  appears  from  the  smoke  pictures  as 
well  as  from  the  results  of  an  analytical  treatment  of  the  wake  of  a  hovering 
rotor  (Reference  4)  that  self-induced  effects  rapidly  carry  the  root  vortices 
up  through  the  center  of  the  rotor  plane  where  they  are  dissipated.  The  tip 
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vortices  are  quite  stable,  however,  and  experience  little  or  no  diffusion  due 
to  viscous  effects  for  several  rotor  revolutions. 

From  the  point  of  view  of  the  fuselage,  the  flow  appears  as  the  superposition 
of  a  steady  free  stream  caused  by  the  translation  of  the  aircraft  and  the 
periodic  flow  induced  by  the  rotor  and  its  wake.  The  contribution  of  the 
fuselage  to  the  flow  at  any  point,  then,  is  essentially  that  due  to  a  body  of 
complicated  geometry  in  unsteady  potential  flow. 


3. 


FORMULATION  OF  THE  ANALYTICAL  MODEL 


If  the  fluid  is  assumed  to  be  inviscid  and  incompressible,  the  flow  described 
in  Chapter  2  has  an  exact  analytical  counterpart.  The  rotor  blades  could  be 
replaced  by  doublet  sheets,  the  wake  by  vortex  sheets  and  the  fuselage  by  a 
surface  distribution  of  sources.  The  strengths  of  these  potential  singularities 
would  be  functions  of  time  and  spacial  coordinates.  The  requirement  of  flow 
tangency  at  the  blades  and  the  fuselage  would  provide  a  pair  of  nonlinear 
integral  equations  for  the  determination  of  the  strength  distributions.  The 
problems  involved  in  the  formulation  and  solution  of  these  equations  preclude 
such  an  approach.  Therefore,  three  major  simplifying  assumptions  have 
been  utilized  in  the  formulation  of  the  model. 


First,  the  rotor  blades  have  been  replaced  by  line  vortices  ox  constant 
strength  spanwise  and  harmonically  varying  strength  azimuthally.  Except 
in  the  immediate  vicinity  of  the  blades,  blade-induced  effects  should  be 
reasonably  well  reproduced  by  such  a  model.  That  is,  more  than  one  or 
two  chord  lengths  from  a  given  blade,  higher  harmonic  variations  in 
circulation  cannot  have  an  appreciable  effect  on  the  flow,  at  least  in  comparison 
with  fuselage  or  wake  effects.  For  present  purposes,  then,  blade  loading  can 
be  regarded  as  a  known  quantity. 


Second,  the  wake  is  represented  by  individual  free  vortices,  one  trailing 
from  the  tip  of  each  blade  vortex;  the  strength  of  each  tip  vortex  varies  along 
its  length  in  accordance  with  the  (known)  azimuthal  variation  of  the  blade 
vortex  from  which  it  trails.  The  wake  vortices  are  assigned  finite  circular 
cores  of  rotational  fluid  (as  physical  vortices  must  have),  so  that  the  velocity 
which  a  vortex  induces  on  itself  may  be  correctly  computed.  The  requirement 
of  constancy  of  circulation  as  specified  by  Kelvin's  theorem  has  not  been  met. 
However,  as  noted  in  Chapter  2,  the  vortex  elements  needed  to  fulfill  that 
requirement  -  shed  vortices  and  trailing  root  vortices  -  do  not  contribute 
significantly  to  the  flow.  For  the  same  reason,  the  detailed  structure  of  the 
wake  near  the  blades,  where  roll-up  occurs,  and  the  vortices  in  the  far  wake 
(more  then  two  or  three  rotor  radii  downstream  of  the  rotor)  have  also  been 
omitted. 
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Lastly,  the  fuselage  is  assumed  to  experience  a  uniform  free  stream  of 
constant  magnitude  and  direction.  This  stream  is  made  up  of  two  components, 
one  being  the  negative  of  the  velocity  of  translation  of  the  aircraft,  and  the 
other  consisting  of  a  time  and  spacial  average  of  the  downwash  which  would 
be  induced  by  the  rotor  and  its  wake  in  the  absence  of  the  fuselage.  At  high 
forward  speed,  the  time  and  spacial  variations  of  the  stream  in  which  the 
fuselage  is  immersed  are  comparatively  small  and  so  may  be  neglected 
without  causing  large  errors.  At  low  forward  speed,  the  downwash  may 
vary  substantially  over  the  fuselage,  but  the  total  effect  of  the  fuselage  on 
the  flow  is  generally  small  compared  to  rotor  and  wake  effects,  so  that  the 
error  is,  presumably,  again  not  appreciable.  It  is  also  assumed  that  the 
flow  about  the  fuselage  is  not  separated,  so  consideration  must  be  limited 
to  regions  where  separation  effects  are  unimportant. 

The  model  for  the  flow  as  so  simplified  can  be  formulated  in  a  straightforward 
manner.  A  coordinate  system  fixed  in  the  tip-path  plane,  as  shown  in 
Figure  l  for  a  two-bladed  rotor,  is  utilized.  The  vector  Vf  in  the  figure 
is  the  free  stream  velocity  (the  negative  of  the  translational  velocity  of  the 
aircraft),  inclined  at  an  angle  aT  to  the  tf-axis  and  parallel  to  the  plane. 
The  angle  V  denotes  the  azimuth  position  of  the  rotor,  and  £  is  the  vector 
positioning  a  given  point  with  respect  to  the  origin. 

The  fluid  velocity  V  at  a  given  point  located  by  the  vector  rp  may  be  expressed 
in  the  form 


i{ip) 


-/  f  F(r)p_f  x  air 

Htt  J  T5 
1 


dA'  +  Vf 


(1) 


where  r,=  rp  -r  ,  F  is  the  local  circulation  about  a  vortex,  and  is  the  source 
strength  per  unit  area  on  the  surface  replacing  the  fuselage.  The  first 
integral  is  a  line  integral  taken  over  all  vortices  in  the  flow,  which  curves 
are  collectively  denoted  by  Cv,  and  the  second  is  an  integral  over  the  surface  AF 
replacing  the  fuselage. 
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The  circulation  r  about  the  blade  vortices  is  known  directly  in  terms  of  flight 
parameters.  The  circulation  about  a  wake  vortex  at  a  given  point  is  simply 
that  value  assigned  to  the  blade  vortex  when  it  was  generating  that  segment 
of  wake.  The  source  strength  <r  may  also  be  regarded  as  known  insofar  as 
the  wake  flow  calculations  are  concerned.  That  distribution  is  assumed  to 
be  such  as  to  produce  the  flow  about  a  body  of  known  geometry  in  a  uniform 
free  stream,  so  <r(r)  may  be  determined  from  a  separate  calculation. 


The  only  information  lacking  for  the  complete  specification  of  the  flow  at  a 
given  instant,  then,  is  the  location  of  the  wake  vortices  at  that  instant.  The 
position  of  a  given  point  on  a  wake  vortex  located  by  the  vector  r  is  the  time 
integral  of  the  velocity  experienced  by  that  fluid  particle: 


(2) 


Thus,  even  after  simplifying,  the  flow  can  only  be  obtained  as  the  solution 
of  the  nonlinear  integral  equation  formed  by  the  substitution  of  Equation  (i) 
into  Equation  (2).  A  direct  analytical  solution  is  not  feasible,  but  the 
problem  is  amenable  to  solution  by  numerical  methods  using  a  high-speed 
digital  computer.  The  manner  in  which  the  formulations  of  Equations  (1) 
and  (2)  were  implemented  for  digital  computation  is  discussed  in  Chapter  4. 
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4. 


IMPLEMENTATION  OF  THE  MODEL 


WAKE  POSITIONING  AND  DISPLACEMENT 

For  purposes  of  numerical  analysis,  a  series  of  points  P1t  P2,  *  •  *,  is 
marked  off  on  each  wake  vortex,  as  shown  schematically  in  Figure  2. 


Figure  2  WAKE  REFERENCE  POINT  IDENTIFICATION 


The  points  are  placed  sufficiently  close  together  so  that  for  purposes  of 
computation  of  wake-induced  velocities,  the  segment  {i )  of  wake  between 
points  P  and  P-  t  may  be  considered  as  a  rectilinear  vortex  having  constant 
circulation  along  its  length.  The  wake  configuration  at  any  instant  is  then 
defined  by  the  locations  of  these  points. 


Computations  are  initiated  by  specifying  that  each  of  the  wake  vortices  lies 
on  a  prescribed  curve.  Generally,  the  curve  is  chosen  to  be  a  skewed  helix 
with  pitch  and  skew  selected  from  simple  momentum  considerations.  The 
numerical  analogues  to  Equations  (1)  and  (2)  are  then  performed  by  first 
summing  the  velocity  contributions  of  all  vortex  and  source  elements  in  the 
flow  at  each  reference  point  (Equation  (1)),  and  then  convecting  each  point 
at  the  computed  velocity  for  a  time  interval  At  (Equation  (2)).  The  time 
interval  At  is  chosen  to  correspond  to  a  small  finite  change  AP in  the 


azimuth  position  of  the  blades: 

At  = 


-TL 


(3) 


where  si  is  rotor  angular  speed.  Once  the  new  coordinates  of  each  reference 
point  are  computed,  the  azimuth  of  the  blade  vortices  is  increased  by  AV 
and  the  velocity  computation  is  performed  again.  As  the  blade  vortices  are 
repositioned,  a  wake  vortex  element  is  added  to  the  flow  at  the  tip  of  each 
blade  vortex,  the  added  vortex  having  a  length  of  approximately  RAP ,  where 
R  is  rotor  radius.  A  corresponding  element  is  dropped  from  the  computations 
at  the  downstream  end  of  each  wake  vortex,  so  the  program  is  not  encumbered 
by  a  wake  of  ever  increasing  size.  The  computations  may  be  continued 
indefinitely  in  this  manner,  but  it  has  been  found  that  a  periodic  wake  is 
eventually  established  in  most  cases,  after  which  there  is,  of  course,,  no 
need  to  perform  further  calculations. 


The  total  number  of  wake  vortices  taken  into  account  and  the  magnitude  of 
Ap  determine  the  accuracy  of  the  flow  representation  at  a  given  point.  It 
has  been  found  that,  for  a  two-bladed  rotor,  a  value  for  Ap  of  thirty  degrees 
is  sufficiently  small  to  define  the  time  variations  of  the  flow.  The  number  of 
wake  elements  considered  depends  on  the  region  of  interest  and  the  forward 
speed.  If  the  flow  beneath  the  rotor  plane  is  desired,  then  the  wake  generated 
by  only  a  few  revolutions  of  the  rotor  is  needed.  It  has  been  found,  for 
example,  that  with  an  advance  ratio,  /z  =  — q;  >  of  about  .  25,  jus'c  two 
revolutions  of  wake  sufficiently  define  the  flow  under  the  rotor  plane.  At 
lower  advance  ratios,  the  free  stream  does  not  clear  the  wake  from  under 
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the  rotor  as  rapidly,  so  more  wake  must  be  retained.  With  jjl  =  ,  15,  four 
revolutions  of  wake  are  needed,  and  for  hovering,  seven  or  eight  revolutions 
must  be  retained  to  compute  the  flow  within  about  0.  5 R  of  the  rotor  plane. 


COMPUTATION  OF  VELOCITY  INDUCED  BY  WAKE  AND  BLADES 


The  velocity  induced  at  an  arbitrary  point  P  by  the  vortices  representing  the 

wake  and  blades  is  simply  the  sum  of  the  effects  of  an  array  of  rectilinear 

vortex  segments.  If  denotes  the  velocity  induced  at  P  by  the  element 

i 

between  points  Pt-  and  ,  it  is  found  from  Equation  (1)  that 


V*. 


=  -----  ( cos  9.  -  cos  9 2  ) 

tj.jrh  ’ 


(4) 


where  /"£  is  the  strength  of  the  element  and  Of ,  and  h  are  defined  in 
Figure  3.  The  velocity  is  directed  normal  to  the  plane  containing  PL ,  p.  f  and 

P  . 


P 


Figure  3  GEOMETRIC  RELATIONSHIPS  DEFINING  THE  FLOW 

INDUCED  BY  A  RECTILINEAR  VORTEX  ELEMENT  "1 


The  necessity  for  specifying  that  the  wake  vortices  have  a  finite  core  size  is 
evident  from  Equation  (4).  if  Pis  coincident  with  either  P;  or  P^f  ,  the 
formula  for  is  indeterminate.  When  computing  wake-induced  flow  at  a 

L 

given  wake  reference  point,  the  effect  of  the  two  wake  elements  adjacent 
to  that  point  must,  therefore,  be  computed  by  taking  core  size  into  account 
(if  a  point  lies  outside  the  core  substance  of  a  vortex,  the  formula  for  induced 
velocity  is  the  same  as  for  a  vortex  with  infinitesimal  core  size,  i.  e. , 
Equation  (4)).  Locally  induced  velocities  depend  directly  on  the  curvature 
of  (the  centerline  of)  the  vortex  at  the  point  in  question.  An  approximation 
to  the  actual  curvature  is  obtained  by  assuming  that  the  two  elements  whose 
effects  at  point  P  are  to  be  computed  lie  on  the  circular  arc  which  passes 
through  points  ,  P-  and  P^1  ,  as  shown  in  Figure  4.  Then  it  is 
assumed  that  the  two  vortex  segments  have  circular  cores  of  rotational 
fluid  of  radii  and  a respectively,  and  that  the  vorticity  within  the  core 
is  tangent  to  the  arc  and  of  a  magnitude  which  varies  linearly  with  distance 
from  the  center  of  the  arc.  The  assumption  of  linear  variation  of  vorticity 
would  seem  plausible,  since  this  is  the  case  for  a  vortex  ring.  The  self- 
induced  velocity  Vs .  at  point  P;  -  i.  e. ,  the  contributions  of  the  elements 
adjacent  to  Pt-  -  may  then  be  written  in  integral  form  by  summing  the  effects 
of  circular -arc  vortex  filaments  of  differential  cross-section.  If  it  is  then 
assumed  that  the  core  radii  are  much  less  than  the  radius  of  curvature  & , 
the  integrand  may  be  simplified  and  the  integral  evaluated.  It  is  then  found 
that 


where  and  are  defined  in  Figure  4.  The  self-induced  velocity  is 

directed  normal  to  the  plane  of  the  arc.  The  core  radius  of  a  given  element 
may  be  assigned  on  a  rational  basis  using  energy  considerations,  as 
discussed  later. 
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CORE 

RADIUS  a-i 


Figure  4  GEOMETRIC  RELATIONSHIPS  DEFINING  THE  SELF-INDUCED 
VELOCITY  AT  WAKE  POINT  PL 


THE  FUSELAGE  REPRESENTATION 

As  noted  previously,  the  source-strength  distribution  is  determined  for  a 
fuselage  in  a  uniform,  steady  free  stream.  Once  this  distribution  is  obtained, 
the  effect  of  the  fuselage  is  taken  into  account  in  a  straightforward  manner 
according  to  the  second  integral  in  Equation  (1). 

The  computation  of  the  source-strength  distribution  may  be  treated  as  a 
separate  problem,  since  the  source  strength  per  unit  area  cr  at  a  given 
pointer  may  be  normalized.  That  is,  cr  may  be  written  in  the  form 

rte)  *  v'*a#cr*(r)  +  cry(r)  +  <rf  (r)  (6) 
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where  v.  ,  V,  and  V.  are  the  components  of  the  free  stream  experienced 

*•<50  VOO  fOO 

by  the  fuselage,  and  the  normalized  strengths  <j^  and  cr^,  are  functions  only 
of  fuselage  geometry.  The  values  for  the  free  stream  velocity  components 
may  be  obtained  by  simply  computing  the  flow  in  the  absence  of  the  fuselage 
and  taking  appropriate  averages. 

The  normalized  source  strengths  are  computed  according  to  the  method 
reported  in  Reference  5.  The  fuselage  surface  is  approximated  by  an  array 
of  quadrilateral  elements.  The  source  strength  per  unit  area  on  a  given 
quadrilateral  is  taken  to  be  constant.  Thus,  if  there  are  Nf  of  these  elements 
replacing  the  fuselage,  there  are  ZNf  strengths  to  be  determined.  The 
requirement  of  flow  tangency  at  Nf  points  (one  on  each  of  the  quadrilaterals) 
then  produces  ZN^  linear  algebraic  equations  for  determining  the  3a£  unknown 

.1. 

strengths.  A  digital  computer  program  was  coded  for  determining  the  source 
strengths  in  this  manner. 

The  method  as  developed  in  Reference  5  was  used  to  obtain  surface  pressure 
distributions,  so  no  assessment  was  made  there  as  to  the  number  of  quadri¬ 
lateral  source  elements  needed  to  adequately  reproduce  the  flow  an  arbitrary 
distance  from  the  fuselage.  Calculations  were  therefore  performed  for  an 
ellipsoid  of  revolution  having  a  ratio  of  minor  to  major  axis  of  0.  4  so  that 
the  numerical  approximation  to  the  flow  could  be  compared  with  an  analytically 
prescribable  flow.  Representation  by  56,  136,  and  184  surface  elements 
were  considered.  The  flows  produced  in  these  three  cases  were  then  compared 
with  the  exact  potential  solution.  The  results  are  shown  plotted  in  Figures 
5  and  6.  Figure  5  compares  the  surface-element  and  exact  solutions  for  the 
free  stream  parallel  to  the  axis  of  revolution,  and  Figure  6  presents  th: 


In  applying  the  method  to  the  flow  about  the  fuselage,  [/  was  assumed 
zero  and  the  fuselage  assumed  to  have  one  plane  of  symmetry,  which 
reduces  the  number  of  unknowns  to  A/f . 
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5  COMPARISON  OF  THE  EXACT  POTENTIAL  SOLUTION  FOR  THE  FLOW  ABOUT 
AN  ELLIPSOID  WITH  APPROXIMATIONS  BY  A  FINITE  NUMBER  OF  SOURCE 
SHEETS  -  AXIAL  FLOW  CASE 
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*  * 
1  results  for  a  cross -flow.  The  quantities  u,  and  plotted  are  the  fluid  velocity- 
components  in  the  %  and  j  directions,  respectively,  for  the  coordinate  system 
|  sketched  on  the  figures. 

t  In  the  case  of  axial  flow,  the  differences  among  the  results  for  56,  136,  and 

*  184  elements  were  too  small  to  be  seen  on  the  plot,  so  only  the  184-element 

case  is  shown  in  Figure  5,  It  may  be  noted  from  this  figure  that  the 
{  approximate  solutions  only  differ  from  the  exact  solution  very  near  the  body, 

and  the  differences  which  do  occur  are  very  small. 

1 

Somewhat  larger  differences  are  apparent  in  the  case  of  cross -flow,  Figure  6, 
j  with  noticeable  improvements  resulting  by  increasing  the  number  of  surface 

i 

i  elements.  The  largest  differences  are  seen  to  occur  along  the  ^.-axis. 

These  are  caused  by  the  singular  nature  of  the  flow  very  near  the  juncture 
.  of  two  surface  elements. 


These  results  would  indicate  that  a  fuselage  with  approximately  the  proportions 
of  the  ellipsoid  considered  would  be  adequately  represented  by  136  surface 
elements.  A  more  accurate  representation  would  not  be  justified  m  the 
light  of  the  errors  introduced  by  assuming  uniform  steady  flow.  It  should 
be  noted,  though,  that  areas  having  large  curvature  must  be  represented  by 
elements  smaller  in  size  than  are  adequate  elsewhere,  thereby  increasing 
the  number  of  elements  required.  Thus,  even  if  the  fuselage  has  approximately 
the  proportions  of  the  ellipsoid  treated,  considerably  more  than  136  elements 
may  be  needed  for  an  accurate  representation  if  the  surface  is  at  all  irregular. 

DETERMINATION  OF  SYSTEM  PARAMETERS 


The  formulations  of  the  model  were  nondimensionalized  for  purposes  of 
coding,  with  lengths  made  dimensionless  by  rotor  radius  R  and  velocities 
by  rotor  tip  speed  SIR.  The  flight  conditions  of  the  aircraft  being  represented 
relate  to  the  computer  program  through  the  following  dimensionless  parameters: 


advance  ratio 


_  Kr 


nondimensionalized  blade  vortex  circulation 


SIR*  > 

tip-path  plane  angle  ocT  and  the  ratio  of  wake-vortex  core  radius  to  rotor 

The  quantities  which  determine ^  are,  of  course,  generally 


radius,  iL 
R 
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available.  Specification  of  the  other  parameters  is  accomplished  in  the 
manner  discussed  below. 


Blade  Vortex  Strength,  P8  (^) 

The  circulation  about  a  rotor  blade  at  a  given  azimuth  position  ^  is  actually 
varying  in  the  radial  direction  (more  or  less  elliptic  ally).  It  was  felJ, 
though,  that  a  sufficient  approximation  to  blade -induced  effects  would  be 
obtained  by  making  the  circulation  constant  spanwise  and  of  a  magnitude 
equal  to  the  maximum  attained  on  the  actual  blade.  It  might  be  supposed 
that  blade-induced  effects  would  then  be  somewhat  exaggerated,  but 
comparisons  of  computed  and  measured  flows  indicate  that  the  error  so 
introduced  is  not  large  enough  to  be  discernable  (see  Chapter  5). 

If  the  circulation  is  assumed  to  vary  elliptically  in  the  spanwise  direction, 
the  maximum  circulation  rm  about  a  blade  is  given  by  (neglecting  root 
effects  and  making  small-angle  approximations  where  appropriate), 

_ 8LbW 

m  rr^ortR2  (  1  sin  ifj)  2 

v/j-aR* 

where  LB  is  the  total  lift  acting  on  the  blade  and^  is  the  air  density.  Then, 
assuming  that  the  total  lift  acting  on  the  rotor  is  constant  and  equal  to  the 
gross  weights  of  the  aircraft,  it  follows  that 

rm 

~  Ztt ?\{1  -  Z sen  </s)  (7) 

where  A  is  a  blade  loading  parameter,  defined  by 

4 IW 

A  "  ttzN8  p  nz  R  + 

with  denoting  the  number  of  blades  on  the  rotor. 
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Inclination  of  the  Tip-Path-Plane,  a.T 


The  tip-path-plane  angle  is  a  complicated  function  of  flight  conditions,  rotor 
dynamics  and  blade  cyclic  and  collective  pitch  settings.  Various  methods 
a^e  available  for  computing  this  angle  (see,  for  example,  Reference  6),  but 
their  use  requires  iterative  or  graphical  procedures  and  the  values  of 
parameters  which  are  generally  not  well  defined. 


Reasonable  values  for  <Xr  have  been  obtained  without  resorting  to  elaborate 
methods  by  using  the  formula 


ccr  = 


2  C[>£ 
TTNq  3 


(8) 


where  C D  is  the  fuselage  drag  coefficient,  defined  as  the  ratio  of  fuselage 
drag  to  3~/pVpZRz  .  Equation  (8)  was  obtained  by  neglecting  rotor  in-plane 
drag  and  fuselage  lift.  The  angle  ocr ,  if  assumed  small,  must  then  be 
simply  the  ratio  of  fuselage  drag  to  rotor  thrust  (aircraft  weight).  Values 
for  ccT  obtained  from  Equation  (8)  were  found  to  be  within  20  to  30  percent 
of  measured  values  for  a  UH-1A  helicopter  (from  Reference  7)  using  a 
value  for  CDf  of  .014  (Reference  8). 


Wake-Vortex  Core  Radius,  a, 


The  nondimensionalized  core  radius  of  each  vortex  element  must  be  specified 
as  it  is  generated  by  a  blade  vortex.  This  core  size  should  correspond  to 
that  of  the  tip  vortex  formed  in  the  physical  flow  when  the  wake  rolls  up. 

In  order  to  obtain  an  estimate  of  the  core  size,  three  major  simplifying 
assumptions  have  been  made.  First,  all  vorticity  in  the  flow  except  that 
contained  in  the  wake  under  analysis  has  been  neglected.  This  may  be 
rigorously  justified  as  an  adjunct  to  the  development  which  follows,  but 
the  calculations  are  somewhat  lengthy  and  so  have  been  omitted.  The 
justification  uses  the  fact  that  some  distance  from  the  wake  the  flow  before 
and  after  roll-up  is  effectively  the  same.  Second,  it  is  assumed  that  the 
wake  may  be  considered  as  though  it  were  generated  by  a  wing  of  the  same 
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geometry  as  the  blade  in  steady  translational  flight.  That  is,  the  tip  and 
root  vortices  of  a  rotor  blade  are  assumed  to  have  the  same  core  size  as 
the  tip  vortices  generated  by  a  wing  of  the  same  aspect  ratio  in  steady 
translational  flight,  and  having  the  same  loading  as  the  blade  had  when  it 
generated  the  segment  of  wake  being  analyzed.  Third,  it  is  assumed  that 
the  roll-up  of  the  wake  of  a  wing  in  steady  translational  flight  can  be  treated 
as  an  unsteady  two-dimensional  problem  rather  than  a  steady  three-dimensional 
one.  The  only  justification  (other  than  intuitive  argument)  for  the  latter  two 
assumptions  is  that  they  do  yield  results  which  appear  to  be  realistic. 

Consider,  then,  a  wing  of  large  aspect  ratio  in  steady  translational  flight, 
as  shown  in  Figure  7.  If  f  denotes  the  circulation  about  the  wing,  then  at 
some  position  x>  in  the  wake  upstream  of  the  point  where  roll-up  begins 
(View  A-A,  Figure  7),  the  circulation  per  unit  length  /(y)  is  given  by 

(9) 

where  ou[y)  is  the  magnitude  of  the  vorticity  in  the  wake  and  My)  is  the 
thickness  of  the  wake  sheet  (see  Figure  7).  Given  the  loading  on  the  wing  - 
i.  e. ,  on  the  rotor  blade  at  some  azimuth  -  y(y)  may  then  be  specified 
through  Equation  (9).  Note  that  /=  0  at  some  point,  say  yQ,  where  r  is  a 
maximum.  Also,  since  r[-b)  =  r(b)  =  0,  -b  corresponding  to  blade  root 
and  b  to  blade  tip,  it  follows  that 

n^O  n  b 

r(<j)cCy  =  J  y(y)ciy  =  rm 
-b  uyQ 

where  fm  denotes  the  maximum  value  of  My).  Replacing  variation  in  x>  in 
the  wake  by  variation  in  time  in  a  plane  parallel  to  the  plane  and  fixed 
in  the  wake,  the  sheet  of  vorticity  would  be  observed  to  roll  up  into  a  pair 
of  vortices  whose  cores  consist  of  the  rotational  fluid  which  was  contained 
in  the  sheet. 

This  process  will  be  assumed,  in  what  follows,  to  be  two-dimensional.  That 
is,  the  flow  in  the  wake-fixed  plane  is  taken  to  be  the  same  as  though  the 
wake  were  extended  to  positive  and  negative  infinity  in  the  ^-direction,  and 
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Figure  7  REPRESENTATION  OF  THE  WAKE  BEFORE 
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that  in  each  plane  parallel  to  that  plane  the  flows  are  identical;  or  taken  still 
another  way,  it  is  assumed  that  derivatives  with  respect  to  x  are  small  in 
comparison  to  changes  in  y  or  insofar  as  roll-up  is  concerned. 

If  the  cores  of  the  vortices  formed  by  the  roll-up  are  assumed  to  be  of 
circular  section  and  of  approximately  uniform  vorticity,  and  if  it  is  further 
assumed  that  co  is  constant  in  the  sheet  before  roll-up  {except  at  ya,  where  it 
changes  sign),  it  follows  that  the  core  radii  of  the  two  vortices  will  be  equal. 
Further,  because  the  moment  of  impulse  of  the  flow  is  conserved,  the 
vortices  locate  on  the  y- axis  with  their  respective  positions  y  =  yL  and 
y  -  yK  given  by 

«<■  ~~  rhzi 

-b 

y*  =  jr  f  *10) 

The  situation  after  roll-up  is  as  shown  in  Figure  8. 


Figure  8  REPRESENTATION  OF  THE  WAKE  AFTER  ROLL-UP 
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If  viscous  effects  are  neglected,  the  kinetic  energy  of  the  flow  per  unit 
depth,  T  ,  must  be  conserved.  By  deriving  expressions  for  T  for  before  and 
after  roll-up  and  equating  them,  a  relationship  for  core  radius  a>  may  be 
obtained. 


From  Reference  9,  T  is  given,  in  general,  by 


where 


T=  - 


P 


>«■» 


Vz 


cLrid£ 


(ID 


The  integrals  are  taken  over  the  entire  area  where  co  is  not  zero.  If  2L  is 

h 

substituted  for  u)  in  Equation  (11)  and  h  is  assumed  to  be  much  less  than  b  , 
the  integrands  may  be  simplified  and  the  integrals  evaluated.  It  is  then 
found  that,  in  the  limit  of  vanishing  sheet  thickness,  T  before  roll-up  is 
given  by 


r=  •  (i2> 


Similarly,  it  is  found  that  after  roll-up, 

T .  ft rA 

2tt 

so  that  the  core  radius  a,  is  determined  from 


+  i 


(13) 


.  l  fb[b 

zu,  J_l  rz 

-b  b  m 


y ~_n 

b 


<Ly  dr)  -  —  . 


(14) 


The  result  of  Equation  (14)  may  be  put  in  a  form  more  directly  related  to 

blade  circulation  as  follows.  Let  r  be  expanded  in  a  Fourier  series, 

according  to  1 

( 

-1 
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(15) 


r(9)  2  ?  ua 

=  —  Z,  G'k  smkG 


r  77" 


where  COsQ  -  .  Then 

b 


rco)  = 


ij^g/7  0  d& 

zrm 


irbscnQ 


E  ka^k  cos  ^  £ 


k--t 


Also,  it  is  found  from  Equations  (10)  and  (15)  that 

■  b 


^  -tkl  *r < 


™  -b 
=  CO, 

If  y  and  (jR  -tjL  are  then  substituted  in  Equation  (14)  and  the  integrations 
performed,  it  is  found  that 


Qj_ 

R 


(16) 


where  r0  is  the  radial  distance  from  the  axis  of  revolution  to  the  blade  root. 

An  indication  of  the  order  of  magnitude  of  core  radius  predicted  by  Equation 

{ 1 6)  can  be  obtained  by  considering  the  case  of  an  elliptic  distribution  for  V 

(i.  e. ,  minimum  induced  drag,  and  hence  maximum  core  radius).  Supposing 

r0  to  be  zero,  and  substituting  =  Ci3  =  •  •  •  =  0,  find  a  value  for 

Bi  of.  0862. 

R 

Some  numerical  computations  were  performed,  utilizing  Equation  (16),  to 
determine  the  order  of  magnitude  of  core  size  predicted  by  this  method  and 
to  provide  an  indication  of  the  amount  by  which  the  core  size  varies  with 
azimuth  and  forward  speed.  Data  for  the  computations  were  generated  by  a 
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blade-loads  program  that  CAL  developed  for  the  U.  S.  Army  Aviation 
Materiel  Laboratories  (these  results  are  generally  in  good  agreement  with 
experimentally  obtained  loadings).  Three  flight  conditions  were  considered: 

(1)  a  UH-1A  rotor  at  105  kt.  forward  speed,  (2)  a  UH-1A  rotor  at  30  kt. 
forward  speed,  and  (3)  a  5 -foot-diameter  experimental  rotor  in  hovering 
flight.  The  ratio  of  rotor  radius  to  vortex  core  radius  was  calculated  for 
four  azimuth  positions  of  the  blades  in  each  case,  using  a  12-point 
harmonic  analysis  of  circulation  distribution.  The  results  of  the  computations 
are  tabulated  below. 


Flight  Condition 

Ratio  of  Rotor  Radius  to  Core  Radius 

ii 

o 

o 

o 

o 

ii 

ii 

i — * 

CO 

o 

o 

V  =  270° 

105  kt.  -  UH-1A 

24.0 

13.  3 

17.  5 

21.  0 

30  kt.  -  UH-1A 

18.  7 

15.  0 

14.  4 

24.  1 

Hover  -  5  Ft.  Rotor 

21.9 

21.  9 

21.  9 

21. 9 

This  method  of  predicting  vortex  core  radius  appears  to  provide  a  reasonable 
order  of  magnitude  for  that  quantity  -  from  about  4.  0  to  7.5  percent  of  rotor 
radius.  Most  notable  is  the  lack  of  any  substantial  variation  in  core  radius 
either  with  forward  speed  or  azimuth  position.  The  variations  obtained  may 
be  considered  small  since  a  change  of  20  or  30  percent  in  core  radius 
produces  a  much  smaller  change  in  self-induced  velocity,  v/ith  the  latter 
varying  as  the  logarithm  of  the  former. 

It  was,  therefore,  possible  to  incorporate  a  major  simplification  in  the 
digital  computer  program,  by  assigning  a  representative  value  of  ,  05  to 
the  ratio  of  core  radius  to  rotor  radius,  independent  of  azimuth  position  or 
flight  conditions.  It  was  verified  th?t  this  simplification  does  not  cause  an 
appreciable  error. 
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5. 


RESULTS  OF  COMPUTATIONS 


GENERAL  REMARKS 


Wake  flow  calculations  were  carried  out  in  two  phases.  The  program  was 
first  operated  with  the  fuselage  representation  omitted.  The  sensitivity  of  the 
wake  model  to  various  parameters  and  the  general  character  of  the  solutions 
were  investigated,  and  a  comparison  of  computed  flow  with  experimental 
measurements  was  made.  Once  it  was  determined  that  the  model  was  both 
operationally  practical  and  consituted  a  valid  representation  of  the  flow, 
extensive  calculations  were  performed  which  incorporated  the  fuselage 
representation.  The  latter  results  correspond  to  a  UH-1B  helicopter  at 
five  different  forward  speeds  ranging  from  hover  to  high-speed  cruise. 

Before  discussing  the  details  of  the  results  obtained,  though,  the  following 
remarks  concerning  the  qualitative  nature  of  the  flows  computed  should  be 
made. 

There  was  some  concern  during  the  formulation  as  to  the  numerical  and/or 
physical  stability  of  the  model;  i„  e.  ,  it  was  questioned  whether  the  approx¬ 
imations  made  in  integration  or  possibly  the  nature  of  vortical  flows  them¬ 
selves  would  allow  a  periodic  flow  to  be  established.  It  was  found,  however, 
that  the  model  generally  does  produce  a  periodic  flow,  even  though  the  initial, 
or  starting,  wake  configuration  may  differ  greatly  from  the  periodically 
varying  wake  geometry.  There  was  some  evidence  of  nonperiodic  behavior 
at  very  low  advance  ratios,  however,  which  appears  to  be  attributable  to  an 
inherent  instability  in  free-vortex  flovs.  At  an  advance  ratio  below  about 
.  07,  wake-induced  effects  are  comparable  with  the  magnitude  of  the  free 
stream,  so  the  individual  vortices  then  interact  to  a  great  extent.  It 
results  that,  for  very  low  advance  ratios,  an  increasingly  nonperiodic  flow 
is  evident  in  proceeding  down  the  wake,  starting  from  about  two  to  three 
revolutions  of  wake  below  the  rotor.  The  wake  geometry  and  the  flow  is 
periodic  above  this  point  (which  is  about  half  a  rotor  radius  below  the  rotor 
plane,  for  moderate  loading).  Halving  the  interval  of  integration  or  adding 
wake  elements  does  not  noticeably  alter  this  situation. 
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The  general  character  of  the  wake  geometries  obtained  is  illustrated  in 
Figure  9.  The  tip  vortices  for  a  two-bladed  rotor  are  shown  drawn  in 
perspective  projection  for  the  case  of  /^  =  .  14,  A  =  .00236,  otT  -  2.  3 
degrees  and  A’/'  =  15  degrees.  A  periodic  flow  was  obtained  for  this  case. 

The  individual  wake  elements  were  projected  as  straight  lines  between  the 
wake  reference  points.  Both  the  initial  wake  configuration  and  the  configuration 
at  y  -  0  after  periodicity  was  established  are  shown  so  that  the  extent  of  the 
wake-induced  distortion  may  be  seen.  The  large  vertical  displacements  two 
to  three  rotor  radii  downstream  (which  are  somewhat  exaggerated  by  the 
perspective)  might  in  part  be  attributable  to  the  size  of  the  integration 
increment  taken,  and  so  are  not  necessarily  representative  of  the  flow  in  the 
far  wake.  The  concentration  of  vortices  at  the  lateral  extremes  of  the  wake, 
gi /ing  the  appearance  of  the  wake  behind  a  conventional  wing,  is  clearly 
evident  in  the  figure. 

COMPARISONS  WITH  EXPERIMENTAL  MEASUREMENTS 

In  order  to  properly  assess  the  model  of  the  rotor,  it  was  decided  to  compare 

computed  and  measured  flows.  Satisfactory  quantitative  data  pertaining  to 

instantaneous  fluid  velocities  could  not  be  located.  A  rather  extensive 

wind  tunnel  program  was  carried  out  by  NACA,  however,  to  measure 

average  velocities  in  the  wake  of  a  two-bladed  teetering  rotor  at  moderate 

to  high  advance  ratios.  The  results  of  the  program  are  presented  in  Reference 

10.  The  data  are  presented  in  Reference  10  as  plots  of  —  as  a  function  of 

_ 

the  three  spacial  coordinates,  where  V  is  the  time  average  of  the  induced 
downwash  (i.  e. ,  with  the  free  stream  component  subtracted  off)  normal  to 
the  tip-path-plane,  and  tg  is  the  value  for  that  quantity  as  predicted  by 
simple  momentum  theory.  For  moderate  to  high  advance  ratio, 

=  JiRCT /(Zfx)  where  CT  is  rotor  thrust  coefficient.  It  was  estimated  by 
the  authors  of  that  report  that  the  error  accumulated  in  measuring  PyVg 
was  less  than  about  +  0.  1  at  the  moderate  advance  ratios. 

A  comparison  was  made  of  the  average  of  computed  velocities  with  these 
measurements  for  one  of  the  test  flight  conditions.  For  the  case  chosen, 
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the  advance  ratio  p,  was  0.  14,  the  thrust  coefficient  CT  was  .  00371,  where 
CT  is  defined  as  the  ratio  of  rotor  thrust  to  Tfps^ ,  and  the  inclination  of 
the  tip-path  plane  <XT  was  2.  3  degrees.  The  comparison  of  computed  and 
measured  results  is  shown  plotted  in  Figure.  10  through  14.  The  data  is 
plotted  for  the  coordinate  system  of  Figure  1,  with  the  coordinates  non- 
dimensionalized  by  rotor  radius. 

Figures  10  through  12  show  the  variation  of  v-/v0  with  y  for  values  of  x  of 
-0.  5,  0  and  0.  5  respectively,  with  each  figure  showing  the  variation  at  both 
-  -0.  07  and  ^  =  -0.  2.  The  computed  and  measured  values  are  seen  to 
agree  quite  well  for  all  three  values  of  x.  It  should  be  noted  that  this 
agreement  is  obtained  even  though  the  instantaneous  computed  values  of  the 
downwash  vary  a  great  deal  from  the  mean.  The  extent  of  the  variation  is 
illustrated  in  Figure  13.  In  this  figure,  the  ratio  of  the  induced  downwash 
vto  is  plotted  against  blade  azimuth  for  x  =  0.  5,  and  y  =  0.  5  for  both 
^  =  -0.  07  and  %  =  -0.  2.  The  computed  and  measured  mean  values  are 
also  indicated  on  the  plots. 

Rather  poor  agreement  was  obtained  downstream  of  the  rotor  in  the  plane. 
This  may  be  seen  in  Figure  14,  where  the  variation  of  v-/ir0  with  is  plotted 
for  x>  -  1. 07  and  ^5  =  2.  07,  with  y  -  0.  The  lack  of  agreement  may  be  due  in 
part  to  the  effects  of  the  structure  supporting  the  test  rotor.  These  effects 
were  not  discussed  by  the  experimenters  in  their  report.  Note  that  at  Z  =  0.  5 
(Figure  12),  the  agreement  was  poor  near  y  =  0,  while  this  is  not  the  case 
at  x>  =  -0.  5  (Figure  10),  where  the  rotor  support  could  have  little  effect  on 
the  flow. 
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Fioure  13  COMPUTED  VARIATION  OVER  ONE  PERIOD  OF  THE  INDUCED  DOWNWASH 
AT  TWO  POINTS  BELOW  THE  ROTOR  PLANE 
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FLOW  CALCULATIONS  FOR  A  UH-1B  HELICOPTER 


Extensive  calculations  were  performed  to  determine  the  flow  in  the  vicinity 
of  a  UH-1B  helicopter.  The  aircraft  was  assumed  to  have  a  gross  weight  of 
5675  pounds  (about  two  thirds  of  maximum  gross  weight)  and  a  rotor  speed 
of  31.4  radians  per  second,  giving  a  value  for  the  blade  loading  parameter, 

A  ,  of  .  00209.  A  total  of  five  cases  were  treated,  corresponding  to  forward 
speeds  from  0  to  110  knots.  The  values  of  the  important  parameters  for 
these  cases  are  tabulated  below. 


Advance  Ratio 

M' 

Tip-Path-Plane 

Angle  oy,  degrees 

Number  of  Revolutions 

of  Wake  per  Blade 

0 

0 

7 

.  0732 

.  655 

5 

.  1465 

2.  62 

4 

.  220 

5.  91 

3 

.  269 

8.  84 

2 

In  subsequent  discussions,  the  above  cases  will  be  identified  by  the 
appropriate  value  of  advance  ratio.  In  all  cases,  A f  =  30  degrees,  and 
of  course,  Ng  =  2.  The  fuselage  was  represented  by  a  total  of  192  source- 
sheet  elements.  The  somewhat  simplified  geometry  used  for  the  surface 
replacing  the  fuselage  is  depicted  in  plan  and  elevation  in  Figure  15.  The 
coordinate  system  in  Figure  15,  and  the  one  used  for  the  plots  of  flows 
following  Figure  15,  is  that  of  Figure  1,  nondimensionalized  by  rotor 
radius. 

Flow  calculations  were  performed  within  the  volume  external  to  the  fuselage 
and  having  the  limits 

-1.1  <  %<  -0.1 
-0.¥-S  <  tj  <  +  O.^-S 

-  0>¥5  5  ^  <  -  0./5 
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In  all  five  cases,  a  periodic  flow  was  established  within  this  volume. 

The  data  are  presented  in  Figures  16  though  55,  at  the  end  of  the  Chapter, 
as  plots  of  the  dimensionless  velocity  components  and  ,  which 

relate  to  their  dimensional  counterparts  by  a  factor  hClR,  The  time 
averages  over  one  period  of  VA,  and  are  also  plotted;  these  averages 
being  denoted  by  l/,,  Vy  and  respectively. 

Caution  should  be  exercised  in  comparing  plots  for  different  advance  ratios. 
Scales  and,  in  some  instances,  the  quantities  plotted,  were  adjusted  for 
each  advance  ratio  to  suit  the  data. 

THE  FLOW  AT  HIGH  ADVANCE  RATIOS 

The  flows  at  advance  ratios  of  .  1465,  .  220,  and  .  269  are  characterized  by 

a  lack  of  unsteadiness  over  most  of  the  region  investigated.  Figures  16,  17 

and  18  show  the  variation  of  the  flow  with  #  at  y  =  0  and  =  -.45  for  these 

three  advance  ratios.  The  variations  for  rotor  azimuths  of  0  and  90 

degrees  and  the  time  averages  for  one  period  are  shown  on  each  figure. 

The  variations  with  *  at  other  azimuth  positions  are  not  significantly 

different  from  the  ones  plotted.  The  flow  is  nearly  constant  with  azimuth 

because  the  fuselage  and  free-stream  contributions  dominate  at  high 

advance  ratios.  This  may  be  seen  from  Figure  19,  which  compares  the 

total  flow  with  the  contributions  Vu.  and  K.  of  the  wake  and  blades  to  Va 

y  o  at  y 

and  J/  ,  for  ^  =  .  220  and  ^  =  0.  Wake  and  blade  contributions  to  V#  were 
too  small  to  be  seen  on  the  scale  used. 

The  dominance  of  the  fuselage  and  free  stream  is  further  illustrated  by 
Figures  20  through  36,  on  which  are  plotted  flow  variations  with  #  and  y  in 
selected  regions  at  the  three  highest  advance  ratios.  In  those  regions  where 
there  were  discernable  flow  variations  with  azimuth,  plots  are  shown  for 
six  azimuth  positions.  Otherwise,  only  the  plots  for  f  =  0  and  ^  =  90  degrees 
are  given.  As  might  be  expected,  somewhat  more  effect  of  the  blades 

is  felt  closer  to  the  rotor  plane  (Figures  26,  31,  and  36,  for  example). 

Spacial  variations  are  still  attributable  primarily  to  the  fuselage,  however. 
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The  assumption  that  the  free  stream  experienced  by  the  fuselage  is  constant 
in  time  is  clearly  justified,  since  time  variations  of  the  flow  in  the  region 
occupied  by  the  fuselage  (  J--.  2)  are  negligible.  It  should  be  noted,  though, 
that  a  similar  statement  cannot  be  made  regarding  the  assumption  that  this 
free  stream  is  also  uniform.  Figure  37  shows  the  spacial  variation  of  the 
downwash  at  f  =  0  in  the  region  where  the  fuselage  is  located,  with  the  flow 
due  to  the  fuselage  omitted  (i.  e.  ,  with  all  source  strengths  zero),  for 
jm  -  .  1465.  The  downwash  can  be  seen  to  vary  by  as  much  as  a  factor  of 
two  in  this  region.  While  undoubtedly  the  flow  is  well  reproduced  some 
distance  from  the  fuselage  by  simply  averaging  these  variations  in  computing 
source  strengths  (the  fuselage  was  found  to  cause  only  minor  distortions 
of  the  wake  at  all  five  advance  ratios)  this  may  very  well  not  be  the  case 
within,  say,  0.  IR  of  the  fuselage.  The  error  introduced  by  taking  an 
average  for  the  downwash  in  computing  fuselage  source  strengths  can  be 
properly  determined  only  by  actually  computing  the  source  strengths  needed 
to  satisfy  the  flow  tangency  condition  for  the  spacially  varying  stream 
depicted  in  Figure  37,  and  then  comparing  the  flow  obtained  using  these 
strengths  with  that  obtained  from  strengths  computed  assuming  the  stream 
to  be  uniform.  These  calculations  were  not  performed  because  of  time 
limitations. 

THE  FLOW  FOR  HOVER  AND  LOW  ADVANCE  RATIO 

There  is  no  large  free  stream  component  of  velocity  to  convect  the  wake 
from  beneath  the  rotor  plane  in  the  cases  of  =  0  and  /o  =  .  0732,  so  the 
opportunity  for  interactions  among  the  wake  elements  is  much  greater  than 
at  high  advance  ratio.  As  a  result,  the  flow  is  correspondingly  altered. 

First,  consider  the  flow  for  =.0732.  The  computational  results  for  this 
advance  ratio  are  plotted  in  Figures  38  through  45.  The  general  character 
of  the  curves  is  seen  to  be  similar  to  that  found  at  the  higher  advance  ratios. 
The  fluctuations  with  azimuth  are  larger,  however.  Close  to  the  rotor 
plane  =  -.  25  and  ^  =  -.  15),  a  net  upwash  is  obtained  in  some  regions,  and 
the  shapes  of  the  curves  are  considerably  altered  with  changes  in  azimuth. 
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The  amount  contributed  to  the  flow  by  the  blades  and  wake  at  /z  -  •  0732  is 
indicated  by  the  plots  of  .Figure  46.  While  the  fuselage  and  free  stream  still 
dominate  the  ^-component  at  this  advance  ratio,  the  flow  normal  to  the  tip- 
path-plane  is  largely  due  to  the  wake. 

The  flow  at  jla.  -  0  differs  greatly  from  that  for  nonzero  advance  ratio  both 
in  the  spacial  variations  and  in  the  magnitudes  of  the  velocity  components. 

The  results  for  hovering  are  plotted  in  Figures  47  through  53.  The  flow 
components  in  the  p  and  u  directions  are  seen  to  be  comparatively  small  - 
on  the  order  of  5  or  10  -  over  the  whole  region  investigated,  while  (4.  varies 
from  small  positive  values  to  negative  values  of  from  30  to  50.  Those  areas 
where  large  downwash  occurs  can  be  readily  identified  as  lying  within  the 
envelope  of  the  wake.  The  large  azimuthal  fluctuations  in  the  flow  are  caused 
by  the  periodic  passage  of  the  wake  vortices  about  this  envelope. 

It  will  be  noted  that  the  downwash  within  the  wake  is  from  three  to  five 
times  greater  in  hover  than  it  is  when  /zis  not  zero,  even  though  the  rotor 
loading  is  the  same  for  all  cases.  An  explanation  for  this  may  be  obtained 
from  simple  momentum  considerations.  The  rotor  thrust  must  be  derived 
from  a  given  net  flux  of  momentum  in  the  negative  ^.-direction,  regardless 
of  the  advance  ratio.  In  hovering  flight,  this  entire  flux  must  be  carried 
by  the  downwash  component  of  velocity.  If  the  aircraft  is  in  forward  flight, 
however,  the  ^-component  of  fluid  momentum  may  be  convected  from  a 
given  control  volume  in  both  the  p  and  ^directions.  Thus,  the  required 
flux  is  obtained  at  nonzero  advance  ratio  without  a  large  downwash  component, 
the  wake  flow  then  being  somewhat  more  comparable  to  that  of  a  wing  than 
a  propeller. 

The  contribution  of  the  fuselage  to  the  total  flow  may  be  seen  from  Figure  54, 
where  the  total  flow  and  that  due  to  the  wake  and  blades  are  compared. 
Variations  with  p  are  shown  for  y  -  .  3,  a.  =  -.  45  and  -  C.  The  fuselage 
may  be  seen  to  alter  the  downwash  by  only  about  ten  percent. 
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In  further  regard  to  the  fuselage,  it  was  found  that  the  use  of  a  uniform,  steady 
downwash  in  computing  the  source  strengths  was  well  justified  for  the  case 
of  /^  =  0.  Figure  55  shows  the  downwash  variation  in  the  z-y.  plane,  at 
2*  =  -.  25  and  ^  =  -»  45  for  y  -  0,  which  results  when  the  fuselage  representation 
is  omitted.  The  downwash  at  other  azimuths  and  over  the  remainder  of  the 
volume  occupied  by  the  fuselage  was  found  to  be  approximately  the  same  as 
shown  in  Figure  55,  i.  e. ,  between  30  and  40.  Thus,  in  hover,  the  fuselage 
lies  well  within  the  wake  envelope  and  is  indeed  subjected  to  a  nearly  constant 
and  uniform  free  stream. 

Also  shown  in  Figure  55  is  the  downwash  in  the  location  of  the  fuselage  for 
juu  -  .  0732.  In  this  case,  large  variations  can  be  seen  to  occur  spacially 
as  wake  cortices  are  swept  directly  over  this  region.  These  variations 
are  much  larger  than  those  which  occurred  at  the  higher  advance  ratios 
(Figure  37).  It  would  appear,  then,  that  the  approximations  used  in 
representing  the  fuselage  are  most  open  to  question  at  the  intermediate 
advance  ratios,  from  about  .  05  to  .  10.  In  this  range,  the  fuselage 
contributes  substantially  to  the  flow  while  experiencing  large  flow  variations 
itself. 
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Figure  26  VARIATION  OF  VELOCITY  COMPONENTS  WITH 
•  FOR  y  =  0,  $  =  -0.  15  AND  jju  -  0.  1465 
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VELOCITY  COMPONENTS  WITH  y  Figure  28  VARIATION  OF  VELOCITY  COMPONENTS 
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Figure  36  VARIATION  OF  VELOCITY  COMPONENTS  WITH  % 
FOR  lj  =  0,  j  =  -0.  15  AND =  0.269 


VARIATION  OF  VELOCITY  COMPONENTS  WITH  x 
FOR  y  -  0.3,  £  =  -0.^5  AND jj.  =  0.0732 
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Figure  46  COMPARISON  OF  TOTAL  FLOW  WITH  WAKE-  AND  BLADE- INDUCED 
FLOW  FOR  y  0.3,  j-  =  -0.45  AND/^=  0.0732 
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Figure  48  VARIATION  OF  VELOCITY  COMPONENTS  WITH  y 
FOR  po  =  -0.25,  =  -0.45  AND ^  =  0 


Figure  49  VARIATION  OF  VELOCITY  COMPONENTS  WITH  * 
FOR  y  =  -0.3,  £  =  -0.45  AND/^=  0 


Figure  50  VARIATION  OF  VELOCITY  COMPONENTS  WITH  * 
FOR  lj  =  0.3,  £  =  -0.45  AND/4  =  0 


Figure  52  VARIATION  OF  VELOCITY  COMPONENTS  WITH  * 
FOR  y  =  0.3,  £  =  -0.25  AND  /x>  -  0 
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Figure  53  VARIATION  OF  VELOCITY  COMPONENTS  WITH  jt> 
FOR  Lj  =  0,  4  =  -0.  15  AND/4=  0 


SPACIAL  VARIATION  OF  DOWNWASH 
IN  THE  ABSENCE  OF  A  FUSELAGE 
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CONCLUSIONS  AND  RECOMMENDATIONS 


On  the  basis  of  the  comparisons  with  experimental  data  and  the  analysis  of 
the  calculations  made  for  a  UH-1B  helicopter,  it  may  be  concluded  that  the 
theory  developed  provides  a  valid  representation  of  the  time-varying  flow  in 
the  vicinity  of  a  helicopter.  The  simplifications  employed  in  formulating 
the  model  exclude  its  application  from  two  specific  regions,  however. 

Because  the  detailed  structure  of  the  wake  in  the  vicinity  of  the  blades  has 
been  omitted,  the  flow  in  the  rotor  plane  itself  cannot  be  accurately 
reproduced.  Similarly,  because  it  was  assumed  that  the  fuselage  experiences 
a  uniform,  constant  free  stream,  the  model  must  be  applied  with  caation  in 
the  immediate  vicinity  of  the  fuselage.  Two  areas  for  refinement  of  the  theory 
then  suggest  themselves. 

In  order  to  allow  for  the  prediction  of  the  flow  in  the  rotor  plane,  some 
approximation  to  the  actual  wake  within  five  or  ten  chord -lengths  of  the 
blades  could  be  added.  For  example,  a  sheet  of  vorticity  of  the  proper 
dimensions  and  strength  distribution  might  be  used.  This  could  then  lead 
to  a  determination  of  blade  loads,  by  computing  the  instantaneous  circulation 
about  the  blade  vortices  from  the  downwash  computed  at  that  instant. 

The  fuselage  representation  could  be  improved  in  two  different  ways.  One 
possibility  is  to  assume  that  the  free  stream  to  which  the  fuselage  is  subjected 
is  constant  in  time  but  varies  spatially,  with  the  variation  determined  by 
computing  the  time  average  of  the  flow  obtained  in  the  absence  of  the  fuselage. 
This  could  be  accomplished  with  only  minor  modifications  to  the  existing 
computer  programs,  and  would  measurably  improve  the  flow  representation 
near  the  fuselage  at  high  advance  ratios. 

The  other  possibility  for  refinement  of  the  fuselage  representation,  while 
requiring  more  extensive  reprogramming,  would  insure  that  the  fuselage 
is  precisely  represented  regardless  of  the  flow  to  which  it  is  subjected. 

This  entails  the  construction  of  an  exact  numerical  analogue  to  the  flow. 
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That  is,  at  each  azimuth  position  of  the  blade  vortices,  the  strengths  of  the 
fuselage  source  elements  would  be  recomputed  so  that  the  flow  at  each 
fuselage  element  is  tangent  to  the  surface;  the-fuselage  boundary  condition 
would  then  be  satisfied  exactly  (within  the  limits  of  the  digital  analogue 
approximation  itself)  both  spacially  and  in  time.  Preliminary  estimates 
indicate  that  this  procedure  would  be  operationally  feasible,  computer 
running  time  only  being  increased  by  a  factor  of  about  two  for  a  representative 
case. 
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